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Abstract. In this paper we develop an axiomatic setup for algorithmic homological 
algebra of AsELian categories. This is done by exhibiting all existential quantifiers en- 
tering the definition of an ABELian category, which for the sake of computability need 
to be turned into constructive ones. We do this explicitly for the often-studied example 
ABELian category of finitely presented modules over a so-called computable ring R, i.e., a 
ring with an explicit algorithm to solve one-sided (in)homogeneous linear systems over R. 
For a finitely generated maximal ideal m in a commutative ring R we show how solving 
(in)honiogeneous linear systems over Rm can be reduced to solving associated systems 
over R. Hence, the computability of R implies that of Rm- As a corollary we obtain the 
computability of the category of finitely presented i?m-modules as an ABELian category, 
without the need of a MORA-like algorithm. The reduction also yields, as a by-product, 
a complexity estimation for the ideal membership problem over local polynomial rings. 
Finally, in the case of localized polynomial rings we demonstrate the computational ad- 
vantage of our homologically motivated alternative approach in comparison to an existing 
implementation of Mora's algorithm. 
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1. Introduction 



As finite dimensional constructions in linear algebra over a field k boil down to solving 
(in) homogeneous linear systems over k, the GAUSSian algorithm makes the whole 
theory perfectly computable, provided k itself is. Solving linear systems in GAUSSian 
form, i.e., in reduced echelon form, is a trivial task. And computing the GAUSSian form of 
a linear system is thus a major step towards its solution. 

Homological algebra of module categories can be viewed as linear algebra over general 
rings. Hence, in analogy to linear algebra over a field one would expect that solving linear 
systems would play an important rule in making the theory computable. [BROS] introduced 
a data structure for additive functors of module categories useful for an efficient computer 
implementation. Solving linear systems was used to describe the calculus of such functors in 
a constructive way. Here we proceed in a more foundational manner. We show that solving 
linear systems is, as expected, the key to the complete computability of the category of 
finitely presented modules, merely viewed as an ABELian category [HS97, Rot09, Wei94]. 
In Section 2 we list all the existential quantifiers entering the definition of an ABELian 
category. Turning all of them into algorithms for any given ABELian category establishes its 
computability. This abstract point of view widens the range of applicability of a computer 
implementation along these lines beyond the context of module categories. 

Section 3 addresses the computability of the ABELian categories of finitely presented 
modules over so-called computable rings. A ring R is called computable if one can 
effectively solve (in)homogeneous linear systems over R (cf. Def. 3.2). Proposition 3.1 
together with Theorem 3.4 show that, as expected, the computability of the ring together 
with some simple matrix operations indeed suffice to provide all the algorithms needed to 
make this category computable as an ABELian category. 

One way to verify the computability of a ring is to find an appropriate substitute for 
the GAUSSian algorithm. Fortunately such substitutes exist for many rings of interest. 
Beside the well-known Hermite normal form algorithm for principal ideal rings with 
computable gcd's, it turns out that appropriate generalizations of the classical GrOBNER 
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basis algorithm for polynomial rings [Buc06] provide the desired substitute for a wide class 
of commutative and noncommutative rings [Lev05, Rob06]. 

Although finding a substitute for the GAUSSian algorithm, which we will refer to as 
computing a "distinguished basis^", is the traditional way to solve linear systems over 
rings, it is only one mean to this end. Indeed, other means do exist: 

Let Ryn be the localization of the commutative ring i? at a finitely generated maximal 
ideal m < R. Theorem 4.1 together with Corollary 4.2 show how computations in the 
ABELian category of finitely presented modules over the local ring Rm can be reduced to 
computations over the global ring R. In particular, one can avoid computing distinguished 
bases over the local ring R^. The idea is very simple. Elements of the local ring R^ can be 
viewed as numerator-denominator pairs (n, d) with r & R and (i G -R\m. Likewise, (r x c)- 
matrices over i?^ can be viewed as numerator-denominator pairs {N, d) with G R^^'^ 
and (i G -R \ m. It is now easy to see that solving (in) homogeneous linear systems over R^ 
can simply be done by solving associated systems over R, thus deducing the computability 
of Rm from that of R. 

In principle, MORA's algorithm, which provides a "distinguished basis", can replace 
BuCHBERGER's algorithm for all sorts of computations over localized polynomial rings. 
This also seems to be the common practice. Nevertheless, a considerable amount of these 
computations only depend on the category of finitely presented modules over localized 
polynomial rings merely being ABELian. From this point of view we show how BuCH- 
BERGER's algorithm suffices to carry out all such constructions and explain in §4.4 why 
our homologically motivated approach to localization of polynomial rings is computation- 
ally superior to an approach based on MORA's algorithm. However, MORA's algorithm 
remains indispensable when it comes to the computation of HiLBERT series of modules over 
such rings, for example. Still, these modules are normally the outcome of huge homological 
computations, which often enough only become feasible through our alternative approach. 
Serre's intersection formula is a typical example of this situation (cf. Example 6.3). 

In Section 5 we will shortly describe our implementation. The examples in Section 6 
illustrate the computational advantage of our alternative approach. An existing performant 
implementation of MORA's algorithm fails for these examples. 

The paper suggests a specification which can be used to realize a constructive setting for 
the homological algebra of further concrete ABELian categories. Realizing this for other 
ABELian categories is work in progress. 

2. Basic Constructions in Abelian Categories 

The aim of this section is to list the basic constructions of an ABELian category with 
enough projectives, which suffice to build all the remaining ones. In case these few basic 
constructions are computable, it follows that all further constructions become computable 
as well. 

In the list we only want to emphasize the existential quantifiers, that need to be turned 
into constructive ones. We decided to suppress the universal properties needed to correctly 



Basis in the sense of a generating set, and not in the sense of a free basis. 
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formulate some of the points below, as we assume that they are well-known to the reader. 
A detailed treatise can be found in Appendix A. 

^ is a category: 

(1) For any object M there exists an identity morphism 1m- 

(2) For any two composable morphisms 0, ip there exists a composition ipip. 

Tl is a category with zero: 

(3) There exists a zero object 0. 

(4) For all objects M,N there exists a zero morphism Omn- 
A is an additive category: 

(5) For all objects M, there exists an addition^ {4',4') ^ (p + ip the ABELian 
group Homyi(M, A^). 

(6) For all objects M, N there exists a subtraction {(pjip) (p ~ ip in the ABELian 
group Homyi(M, A^). 

(7) For all objects Ai, A2 there exists a direct sum Ai © A2. 

(8) For all pairs of morphisms (pi : Ai ^ M, i = 1,2 there exists a coproduct mor- 
phism (01, 02) : Ai © M. 

(9) For all pairs of morphisms (pi : M ^ Ai, i = 1,2 there exists a product morphism 
{01,02} : M ^ Ai © A2. 

A is an ABELian category: 

(10) For any morphism (p : M ^ N there exists a kernel ker0 M. 

(11) For any morphism t : L ^ M and any monomorphism n : K ^ M with r0 = 
for = coker k there exists a Uft tq : L ^ K oi t along k. 

(12) For any morphism (p : M N there exists a cokernel A^ coker 0. 

(13) For any morphism rj : N L and any epimorphism e : N C with (prj = for 
= ker e there exists a colift tiq : C ^ L of rj along e. 

A has enough projectives: 

(14) For each morphism (p : P ^ N and each morphism e : M ^ N with im0 < ime 
there exists a projective lift 0i : P — > M of a along e. 

(15) For each object M there exists a projective hull u : P ^ M. 

If for an ABELian category A we succeed in making the above basic constructions com- 
putable, all further constructions which only depend on A being ABELian will be com- 
putable as well (see for example Appendix B). 

Definition 2.1. Let A be an ABELian category. 

(1) We say that A is computable as an ABELian category if the existential quan- 
tifiers in (1)-(13) can be turned into constructive ones. 

^The last point of Remark A. 6 states that the addition can be recovered from the product, coproduct, 
composition, and identity morphisms. A direct description of the addition is nevertheless important for 
computational efficiency. 
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(2) If additionally the existential quantifiers in (14)-(15) can be turned into constructive 
ones, then we say that A is computable as an ABELian category with enough 
projectives. 

[Bar] details a construction of spectral sequences (of filtered complexes) only using the 
axioms of an ABELian category as detailed above. In particular, all arguments are based 
on operations on morphisms rather than chasing single elements. 

To compute the derived functors of an additive functor F : 71 — )■ !B where A does not 
have enough projectives resp. injectives one needs to provide a substitute for projective 
resp. injective resolutions. The abstract DE Rham theorem suggests the use of so-called 
left (resp. right) F-acyclic resolutions (as used in [Har77, Prop. III. 6. 5], for example). 

3. COMPUTABILITY IN AbELIAN CATEGORIES OF FINITELY PRESENTED MODULES 

The previous section suggests a short path to ensure computability in an ABELian cat- 
egory. This section follows that path for the often-studied example of module categories. 

From now on let ^ := i? — fpmod be the category of finitely presented lejt i?-modules. 
The category fpmod — R of finitely presented right i?-modules is treated analogously. In 
this section we will show how to make the basic operations of §2 computable. As customary 
from linear algebra the basic data structure for computations will be finite dimensional 
matrices over R. 

Finitely presented i?-modules are in particular finitely generated. Thus, a morphism in 
i?— fpmod can be represented by a finite dimensional matrix, the so-called representation 
matrix, with entries in R. 

A free object'^ in R — fpmod is a free module of finite rank r, i.e., a module of the 
form R^^"^ . And since every finitely generated module M is an epimorphic image of some 
V : Fq ^ M with Fq = _R^^^°, it follows that R — fpmod even has enough free objects. 

By definition of i?— fpmod each object even admits an exact sequence Fi ^Fq — >^M, 

Fi = R^^^^ being free of finite rank. The morphism di is called a finite free presentation 
of M. If we denote by M G R^^^'^o ^]-^g matrix representing di and call it presentation 
matrix of M, then u induces an isomorphism from 

cokerM := i^^^'^Vi^^^'^M = coker(i?i^'^^ 4 i?^^"«) = coker^i 

to M. The rows of M are regarded as relations among the tq generators of M given by the 
residue classes of the unit row vectors in R^^'^'o / R^^'^'^n (cf. [BROS, §2], [GP08, Def. 2.1.23], 
[DL06, Def. 1.11]). 

Denote hy R — fpres the category of finite left /^-presentations with objects being finite 
dimensional matrices over R, where one identifies two matrices M G i?^i^^o and M' G R^i^'^o 
with the same number of columns to one object, if R^^^^}'i = R^^^^Vl', as i?-submodules 
of R^^^°. The set IIomij_fpres (M, L) of morphisms between two objects M G R^'^^'^o and 
L G R'^'^^^'o is the set of all Tq x rg-matrices Lp over R with R^^'^^Kip < R^^^'^L, where one 



^Cf. [HS97, §11.10]. 
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identifies two matrices ipi and (p2 to one morphism, if tliey induce tlie same i?-module 
liomomorpliism from coker M to coker L. Summing up: 

Proposition 3.1. R — fpres ^°'^'^> R — fpmod is an equivalence of categories. 

Tlie advantage oi R — fpres is tliat it can directly be realized on a computer. Hence, 
describing the basic constructions of §2 in i? — fpres makes the category R — fpmod 
computable. 

But we note that for R — fpmod or equivalently R — fpres to be AsELian, cokernels 
and kernels of morphisms between finitely presented modules need to be finitely presented. 
This is obvious for cokernels but in general false for kernels: 

Assumption (*): 

From now on we assume that i? is a ring for which the category R — fpmod is 
AsELian. 

Left NOETHERian rings are the most prominent rings satisfying this assumption. 

3.1. Basic operations for matrices and computable rings. Let A be an ri XTg-matrix 
over R. 

3.1.1. An X G K"^^^^ is called a matrix of generating syzygies (of the rows) of A if 

for all X G R^^^^ with xA = 0, there exists a y G R^^''"^ such that yX = x. The rows of X 
are thus a generating set of the kernel of the map 

X = SyzygiesGenerators(A) 

and say that X is the most general solution of the homogeneous linear system XA = 0. 

Further let L be an r[ x ro-matrix over i?. X G R^"^^"^^ is called a matrix of relative 
generating syzygies (of the rows) of A modulo L if the rows of X form a generating 

set of the kernel of the map R^^''^ A coker L. We write 

X = RelativeSyzygiesGenerators(A, L) 

and say that X is the most general solution of the homogeneous linear system XA + YL = 
0. This last system is of course equivalent^ to solving the homogeneous linear system 




= (cf. [BROS, §3.2]). 



3.1.2. Further let B be an r2 x rg-matrix over R. Deciding the solvability and solving 
the inhomogeneous linear system XA = B is equivalent to the construction of matrices N, T 
such that N = TA + B satisfying the following condition: If the i-th row of B is a linear 
combination of the rows of A, then the i-th row of N is zero^. Hence the inhomogeneous 

^In practice however, one can often implement efficient algorithms to compute X without explicitly 
computing Y. 

^So we do not require a "normal form", but only a mechanism to decide if a row is zero modulo some 
relations. 
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In case N = we write X = RightDivide(B, A). 

Rows of the matrices A and B can be considered as elements of the free module R^^^°. 
Deciding the solvability of the inhomogeneous linear system XA = B for a single row matrix B 
is thus nothing but the submodule membership problem for the submodule generated 
by the rows of the matrix A. Finding a particular solution X (in case one exists) solves the 
submodule membership problem effectively. 

As with relative syzygies we also consider a relative version. In case the inhomogeneous 
system XA = B mod L is solvable, we denote a particular solution by 



Definition 3.2. A ring R is called left (resp. right) computable if any finite dimen- 
sional inhomogeneous linear system XA = B (resp. AX = B) over R is effectively solvable 
in the following sense: There exists algorithms computing SyzygiesGenerators(A) and 
DecideZeroEf f ectively(A, B). i? is called computable if it is left and right computable. 

In other words, a ring R is computable if one can effectively solve (in)homogeneous linear 
systems over R. 

Remark 3.3. If the ring R is left computable then the categories R — fpres and (hence) 
R — fpmod are AsELian, and the Assumption (*) is satisfied. 

We want to emphasize that all the free modules used in the constructions below are 
assumed to be given on a free set of generators. This is necessary since there is no known 
algorithm to decide whether a finitely presented module over a computable ring R is free 
or not. In practice this means that if we need to construct a free left i?- module of rank r 
we simply present it by the empty matrix in i?"^*". 

3.2. Computability of the category of finite presentations. 

Theorem 3.4. Let R be a (left) computable ring and A := i?— fpres the ABELian category 
of finite left R-presentations. Then A is computable as an ABELmn category with enough 
projectives. 

Proof. Using all of the vocabulary introduced so far we show how for the category i?— fpres 
the 15 operations for ABELian categories listed in §2 can be turned into algorithms. 

In the following we denote M G R^^^'^o g^j-^fj ]\j ^ j^sixso presentation matrices of M and 
N, respectively. I.e., M := cokerM and N := cokerN. 

R — fpres is a category: 

(1) IdentityMatrix: The identity morphism 1m of M := coker(i?^^''i A R^^''°) is 
represented by the identity matrix I^o G R^°^^°. 



(M, T) = DecideZeroEf fectively(A, B) and N = DecideZero(A, B). 



X = RightDivide(B, A, L). 



This is equivalent to solving (X 
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matrix 



(2) Compose: The composition of two composable morphisms 0, if) represented by the 
matrices A, B is represented by the matrix product AB. 

R — fpres is a category with zero: 

(3) A zero object is presented by an empty matrix in R^^^ . 

(4) ZeroMatrix: The zero morphism Omtv for pairs of objects M := coker^R — )■ 
^ixro^^ jY ._ coker(i?"^^*^ A R^^^°) is represented by the zero matrix Oro,so G 

R — fpres is an additive category: 

(5) AddMat: The addition of two morphisms 0, : M — )■ iV represented by the matrices 
A, B is represented by the matrix sum A + B. 

(6) SubMat: The difference of two morphisms (p^ip : M N represented by the matri- 
ces A, B is represented by the matrix subtraction A — B. 

(7) DiagMat: The direct sum of two objects M and is presented by the bfock diagonal 
0^ 

(8) UnionOfRows: The coproduct morphism {4>,ip) of two morphisms (f) : M ^ L and 
ijj : N L represented by the matrices A and B is represented by the stacked matrix 

^A^ 

(9) UnionOf Columns: The product morphism {0, ip} of two morphisms (p : L ^ M and 
ip : L ^ N represented by the matrices A and B is represented by the augmented 
matrix (A B). 

R — fpres is an ABELian category: 

(10) (Relative) SyzygiesGenerators: To compute the kernel ker0 A- M of a mor- 
phism (j) : M ^ N represented by a matrix A we do the following: First compute 
X = RelativeSyzygiesGenerators(A, N), the matrix representing k. Then ker0 is 
presented by the matrix K = SyzygiesGenerators (X). 

(11) DecideZeroEf f ectively: Let r : L — )■ M be a morphism represented by a matrix 
B and k : K "-^ M a monomorphism represented by a matrix A with r0 = for 
= coker/t. Then the matrix X = RightDivide(B, A, M) is a representation matrix 
for Tq : L —> K, the lift of r along k. It is an easy exercise to see that X indeed 
represents a morphism (cf. [BROS, 3.1.1, case (2)]). 

(12) UnionOfRows & IdentityMatrix: The cokernel module coker0 of a morphism 
: M — )■ = coker N represented by the matrix A is presented by the stacked 

/a\ 

matrix ( 1 • The natural epimorphism A^ coker is represented by the identity 
matrix I^o e i^^o^"*". 

(13) Without loss of generality assume that the cokernel module C = coker is presented 
according to (12), with I^^ the representation of the natural epimorphism e : A^ -» 
C. Further let 77 : A^ — t- L be a morphism represented by B. Then the colift 
?7o : C — !■ L along e is again given by the matrix B. 
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R — fpres has enough free objects: 

(14) DecideZeroEf f ectively: Let F be a free i?-module presented by an empty matrix, 
i.e., F is given on a set of free generators. Further let (p : F ^ N and e : M — be 
morphisms represented by the matrices B and A, respectively. The image condition 
im0 < ime guarantees the existence of the matrix X = RightDivide(B, A, N), which 
is a representation matrix of a free lift 01 : F — ?► M along e (cf. [BROS, 3.1.1, case 

(I)])- 

(15) IdentityMatrix: A free hull u : F M of M is given hj F = cokerF with 
F = G R^^''^' and u is represented by the identity matrix I^^. 

□ 

In the constructive setting we seek, it is necessary to decide if a module is zero and if two 
morphisms are equal. To decide if M = coker M is zero check whether DecideZero(l,f.(,, M) = 
0. To decide the equality of two morphisms 0, ■?/' : M — )■ = coker N represented by A and 
B check whether DecideZero(A — B,N) = 0. This in turn enables deciding properties like 
monic, epic, exactness of two composable morphisms, etc. 

As i?-modules are sets in makes sense to compute in M = coker M. This in turn requires 
deciding equality of two elements in m,m' G M, represented by two rows m, m' G R^^^°, 
respectively. This is again achieved by test whether DecideZero(m — m',M) = 0. 

3.3. Closed symmetric monoidal ABELian categories. The category of /^-modules 
over a commutative ring R admits further constructions. The tensor product M ®r A^ of 
two i?-modules M, N turns the R — fpmod into a symmetric monoidal category. It is even 
a closed symmetric monoidal category with the homomorphism module Homij(M, A^) as 
an internal Hom object. 

For the constructibility of the tensor product and its derived functors Torf in R — fpres 
see [GP08, Example. 7.1.5] or [DL06, Problem 4.7]. For the (internal) Hom module and 
the higher extension functors Ext^ see [GP08, Example 2.1.26], [DL06, Problems 4.5,4.6], 
[KROO, Thm. 3.3.15]. The morphism part of these and other functors is systematically 
dealt with in [BROS]. 

3.4. Free and projective modules. Whereas a finitely presented module is free on free 
generators if an only if its presentation matrix is zero, deciding the freeness of a finitely 
presented module over a computable ring, let alone computing a free basis, can in general 
be highly non-trivial. Deciding project iveness is often easier than deciding freeness: Let 

: Fq ^ M be a free presentation of the i?-module M. It follows that M is projective if and 
only if 1/ admits a section a : M Fq (i.e., au = id-M)- Finding the section a for a finitely 

and freely presented module M <^ Fq ^ Fi leads to solving a two-sided inhomogeneous 
linear system^ over R, which can of course be brought to a one-sided inhomogeneous linear 
system if R is commutative. Hence, testing projectiveness of finitely presented modules 
over commutative computable rings is constructive [ZL02]. Another simple exercise is to 
see that an i?-module is projective if and only if Ext)j(M, i^'i(M)) = 0, where Ki(M) is 



'X + YM = Id, MX = 0, where X is a square matrix representing a and Y another unknown matrix. 
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the first syzygy module of M. Tliis can also be turned into an algorithm for commutative 
computable rings as outlined in §3.3. Without the commutativity assumption Serre's 
Remark [Ser55] states that a module admitting a finite free resolution is projective if and 
only if it is stably free. Shortening the finite free resolution [Lam99, Prop. 5.11] yields a 
simple proof that can be turned into an algorithm for computing the projective dimension 
of M (and hence to decide project iveness) whenever M is a left (resp. right) module over 
the left (resp. right) computable ring R and a finite free resolution of M is constructihle 
[QR07, Algorithm 1]. For a not necessarily commutative ring R with finite global dimension 
there is yet another approach based on AuSLANDER's degree of torsion-freeness [AB69], 
which involves the higher extension modules with values in R of the so-called AuSLANDER 
dual module of M. This approach is constructive if R is (left and right) computable and 
the finite global dimension of R is explicitly known (cf. [CQR05, Thm. 7]). 

The QuiLLEN-SuSLIN Theorem states that a polynomial ring . . . over a prin- 
cipal ideal domain k is Hermite, i.e., every stably free . . . , x„]-module is free. Algo- 
rithms were given in [LS92, LWOO, GV02] and implementations in [FQ07] and [CQR07]. 
A constructive version of Stafford's Theorem [Sta78] offers a way to decide freeness and 
compute a free basis of finitely presented stably free modules over the Weyl algebra An{k) 
and the rational Weyl algebra Bn{k), where is a computable field of characteristic 
[QR07]. The Jacobson normal form [Coh85] offers an alternative algorithm for Bi[k). 
An implementation can be found in the Maple package Janet [BCG''"03]. 

Computability stands for deciding zero and computing syzygies. As mentioned in the 
Introduction, the axiomatic approach pursued so far underlines the conceptual importance 
of solving (in)homogeneous linear systems rather than computing a "distinguished basis", 
the latter being the traditional way to solve such systems. 

4. Computing over local commutative rings 

This section provides a simple alternative approach to solve (in) homogeneous linear sys- 
tems over localizations of commutative computable rings at maximal ideals. The simplicity 
lies in avoiding the computation of distinguished bases over these local rings. In particu- 
lar, for localized polynomial rings this approach offers a way to circumvent MORA's basis 
algorithm. More precisely: 

For a commutative computable ring R with a finitely generated maximal ideal m we show 
in Lemma 4.3 and Proposition 4.5 how to reduce solving linear systems over the local ring 
Rm to solving linear systems over the "global" ring R. Thus, the computability of R implies 
that of and as a corollary, the computability of R^ — fpres as an ABELian category. 

This reduction leads to a result in complexity analysis of local polynomial rings, formu- 
lated in §4.3. Finally, in §4.4 we compare our approach to MORA's algorithm in the case 
of localized polynomial rings. 

Let 1 G S* C be a multiplicatively closed subset of the commutative ring R. Recall, 
the localization of at S' is defined by S~^R := \ r G R,s E 5*}/^, where ^ ~ 7 ^ 
G {rs' — sr') ■ S. The localization S~^M of an i?-module M is defined by S^^R^n M and 
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the localization S^^(f) of a morphism (p maps y G S~^M to ^ipA, Note that localization is 
an exact functor. For a prime ideal p < R the set S* := i? \ p is multiplicatively closed 
and the localization Mp at p is defined as S~^M. In this section we will only treat the case 
when the prime ideal p = m is maximal. 

4.1. Computability in the category of finite presentations over localized rings. 

Theorem 4.1. Let R be a commutative computable ring and m = (mi, . . . , ruk) a finitely 
generated maximal ideal in R. Then Rm is a computable ring. 

The proof of this theorem will be given in §4.2. We first draw the following conclusion: 

Corollary 4.2. Let R and m as in Theorem 4-1- Then R^ — fpres is ABELian and is 
computable as an ABELzan category with enough projectives. 

Proof. Objects and morphisms in R^^ — fpres are given by finite matrices with entries in 
the localized ring R^, being fractions of elements of R. Each such matrix can be rewritten 
as a fraction * with a numerator matrix A over R and a single element s G S* := i? \ m 
as a common denominator^. Using this simplified data structure for matrices over R^ 
is mandatory for the proof of Theorem 4.1. The computability of R^, essential for the 
constructions (10), (11), and (14), is the statement of Theorem 4.1. We still need to go 
through the remaining basic constructions listed in the proof of Theorem 3.4 and adapt 
the required matrix operations to the new data structure: 

Points (l)-(4) are covered by taking the identity matrix (with 1 as denominator), com- 
posing the matrices ^ and | to ^, taking an empty matrix, and taking the zero matrix 
(with 1 as denominator), respectively. Points (5)-(9), stemming from the axioms of an 
additive category, are also easily seen to reduce to the corresponding constructions over 
the global ring R after writing the involved pairs of matrices with common denominators. 
The presentation matrix and representation matrix of the corresponding natural epimor- 
phism of the cokernel (12) are again given by stacking matrices after writing them with 
a common denominator and by taking the identity matrix, respectively. The colift (13) 
as described above was trivial anyway. The presentation matrix of a free hull (15) of a 
module M generated by tq elements is again given by an empty matrix with tq columns, 
while the identity matrix still represents the natural epimorphism. □ 

4.2. Proof of Theorem 4.1. Let i? be a commutative computable ring. Further let 
m = {nil, ■ ■ ■ ,^k) a finitely generated maximal ideal in R. For the computability of R^ 
we need to compute generating sets of syzygies and to (effectively) solve the submodule 
membership problem (cf. Def. 3.2). This is the content of Lemma 4.3 and Proposition 4.5. 

4.2.1. (Relative)SyzygiesGenerators . A matrix of generating syzygies over R is also a 
matrix of generating syzygies over Rm'- 

'''This leads to a simple data structure for matrices over R^, which will prove advantageous from the 
standpoint of a computer implementation (cf. §5.1). 
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Lemma 4.3 (Syzygies). Let A G i?™""". Rewrite A = | with A G i?™''" an(i S G i? \ m. // 

X G R^^^ is a matrix of generating syzygies for k, then the matrix X := j zs a matrix of 
generating syzygies of k. 

Proof. Starting from an exact sequence — ?■ i?^^^ A R^^"^ \ i?^^" exactness of the local- 
ization functor yields an exact sequence Since multiplication 
with 4 is an isomorphism, the sequence -> Rl""" 4 i?^^"^ A i?^^" is also exact. □ 

This lemma is valid more generally for any multiplicatively closed set 5*. Of course, the 
result of the algorithm can have redundant generators and thus might be a non-minimal 
set of generating syzygies. 

Remark 4.4. Nakayama's Lemma implies that a finitely presented i?m-Hiodule M = 
cokerM is given on a minimal set of generators, if and only if its presentation matrix 
M is unit-free (cf. [CLO05, Chap. 5, Prop. 4.3]). An element ^ G -Rm is a unit iff the 
numerator d is not contained in m, which is an ideal membership problem in R. Some 
elementary matrix transformations can now be used to construct unit-free presentation 
matrices (see for example [BROS, §2]). Another easy consequence of Nakayama's Lemma 
is the characterization of minimal resolutions over Rm as the unit-free ones [Eis95, §19.1]. 
Being able to detect units in the computed syzygies we can use a standard procedure^ to 
compute a unit-free and hence minimal resolutions over 

4.2.2. DecideZero(Ef f ectively) . The algorithm for effectively deciding zero is a bit more 
involved. As seen in the proof of Lemma 4.3, denominators of matrices can be omitted 
since multiplication with them is an isomorphism. So without loss of generality we assume 



all denominators 1. To ease the notation let m := 




denote the column of 



generators of the maximal ideal m. 

Proposition 4.5 (Submodule membership). Let k = j E i?™^" and b = j G Rl^"' with 
numerator matrices A and b over R. There exists a row matrix t G R}^"^ with tA + b = 
iff there exists a matrix s G satisfying 

Proof Write s = (y z) with y G i?^^"^, z G R^^''. 

If s and thus y and z exist, we simply set t := zm + 1 ^ m and t := |. For the converse 
implication write t = | with t G and t E R\m. Since t ^ m it has an inverse y 



'Cf. [BROS, §3.2.1], for example. 
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modulo m, i.e., there exist zi, . . . , & R with yt = ziuii + . . . + ZhTUk + 1. We conclude 

tA + B = tA + tb = 

■x^ ytk + {zirrii + . . . + Zkirik) b + 1 ■ b = 
yA + zmb + b = 

for y := yt and z := (zi ... z^) G R^^^ . □ 

This proof is constructive. Another short but nonconstructive proof was suggested by our 
colleague Florian Eisele using Nakayama's Lemma, which is also the most illuminating 
way to interpret Formula (1). The proof given here is indeed a reduction to the effective 
submodule membership problem over R. We don't see a way to generalize this proof to 
non-maximal prime ideals. 

Note that we formulate the proposition for a single row matrix b. For a multi-row 
matrix B simply stack the results of the proposition applied to each row b of B. Contrary 
to GrOBNER basis methods, the proof of the proposition does not provide a normal form^. 
Nevertheless the proof yields an effective solution of the submodule membership problem 
(see §3.1.2). Further note that at no step we need to compute any kind distinguished basis 
over the local ring. The advantages but also the drawbacks of avoiding such a basis will 
be discussed in §4.4. 

Example 4.6. Let R = k[x] for an arbitrary field k and m = (x). We want to compute 
t G R^^ with tA + b = for A = A := x — x"^ and b = b := x regarded as 1 x 1-matrices. 
In this case the gcd of A and mb coincides with b and the extended EuCLIDian algorithm^^ 
directly yields the desired coefficients y and z: 

( X ■ X ) = -X <^ ( -1) ■ (x - x'') + (-X + 1). X . = 

m b 





4.3. Complexity estimation for local polynomials rings. As E. Mayr has shown in 
[May89], the ideal membership problem over the polynomial ring Q[xi, . . . , x„] is exponen- 
tial space complete. Proposition 4.5 implies a result about the complexity of computations 
over the local polynomial ring Q[xi, . . . ,Xn]{xi,...,x„)'- K gives a polynomial time reduction 
from the ideal membership problem over the localized polynomial ring Q[xi, . . . , Xn]{xi,...,xn) 
to the ideal membership problem over the polynomial ring Q[xi, . . . , 

Corollary 4.7. The ideal membership problem over Q[xi, . . . ,Xn]{xi,...,x„) is solvable in 
exponential space. 



^This is an instance where the name DecideZero makes more sense than Reduce or NormalForm. 
"'^'^This is a special case of the Hermite normal form algorithm applied to a one-column matrix. 
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4.4. A comparison with Mora's algorithm. For a polynomial ring R := k[xi, . . . , Xn] 

over a computable field k Mora's algorithm [Mor82] makes the category — fpmod for 
m = (xi, . . . , Xn) computable. The algorithms suggested in Lemma 4.3 and Proposition 4.5 
describe our alternative approach to establish the computability of R^, whereas MORA's 
algorithm can be seen as the classical way to this end. 

Except for the use of a a different reduction method capable of dealing with a local 
term ordering^^, MORA's algorithm proceeds exactly in the same way as BuCHBERGER's 
algorithm, where of course leading terms and s-polynomial depend on the chosen term or- 
dering. And since MORA's algorithm computes, as a by product, the leading ideal/module 
one can read off HiLBERT series. MORA's different reduction method minimizes the so- 
called "ecart" (which measures the distance of a polynomial from being homogeneous) by a 
kind of elimination procedure, which can be very expensive. This makes MORA's reduction 
slower than BuCHBERGER's. For a modern treatment of the theory of local standard bases 
in polynomial rings we refer to [GP08, §1.6, 1.7]. A free implementation can be found in 
Singular [GPS09]. 

Hence, we argue that our approach to localization which only requires an implementation 
of BuCHBERGER's algorithm is computationally superior to a comparable implementation 
of Mora's algorithm. 

Mora's slower reduction is even dwarfed by yet another major issue: Unlike for small 
input, we experienced that Mora's algorithm does not scale well enough when applied to 
large matrices. In comparison with the reduction given by Proposition 4.5, we observed 
that Mora's algorithm creates larger units in R^, which in subsequent computations have 
to be interpreted as denominators. And when writing several matrices over R^ with a 
common denominator, these large units blow up the entries of the numerator matrices. 
Clearly, this makes succeeding computations much harder. 

Remark 4.8. Mora's algorithm is still indispensable when it comes to computing HiL- 
BERT series (and related invariants), which cannot be computed using our approach. Inter- 
mediate homological computations tend to become huge, even if the final result is typically 
much smaller. Our approach is thus suited to get through the intermediate steps, while 
Mora's algorithm can then be applied to the smaller result, e.g. to obtain invariants. This 
is demonstrated in Example 6.3. 

5. Implementation and Data Structures 

As mentioned in the Introduction, this paper suggests a specification for implementing 
homological algebra of ABELian categories (cf. [CS07]). This specification is realized in 
the homalg project [hpalO]. 

We found the programming language of GAP4 ideally suited to realize this specifica- 
tion. GAP4 provides object oriented and (to some useful extent) functional programming 
paradigms, classical method selection, multi-dispatching, and last but not least so-called 
immediate and true-methods, which are extensively used to teach GAP4 how to avoid un- 
necessary computations by applying mathematical reasoning. All these capabilities build 



Local term ordering implies Xi < 1 for all i. 
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upon a type-system, which is as simple as possible and as sophisticated as needed for the 
purposes of high level computer algebra [BL98]. 

The abstract setting of ABELian categories is implemented in the homalg package [BarlO] 
according to §2. Only building upon the basic constructions as abstract operations, the 
implementation provides routines to compute (co)homology, derived functors, long exact 
sequences [BROS], Cartan-Eilenberg resolutions, hyper-derived functors, spectral se- 
quences (of bicomplexes) and the filtration they induce on (co)homology [Bar], etc. 

In order to use these routines for performing computations within a concrete ABELian 
category, the latter only needs to provide its specific implementation of the basic construc- 
tions. The specifics of our implementation for the category i? — fpmod of finitely presented 
modules over a computable ring R were detailed in §3, utilizing the natural equivalence 
coker : R — fpres ^ R — fpmod. More precisely, the proof of Theorem 3.4 shows how this 
equivalence of categories is used to translate all constructions in the ABELian category 
R — fpmod to operations on matrices. Note that this translation is independent of the 
computable ring i?, a point reflected in our implementation. 

This allows the matrices over speciflc computable rings with all their operations to reside 
outside GAP4, preferably in a system that has performant implementations of all the matrix 
operations mentioned in the proof of Theorem 3.4. In turned out that GAP4 does not need 
to know the content but only few characteristic information about the matrices created 
during the computations, minimizing the communication between GAP4 and the external 
system drastically. For further details the interested reader is referred to the documentation 
of the homalg project [hpalO]. 

5.1. LocalizeRingForHomalg. The algorithms presented in §4 are implemented in a GAP4- 
package [GAP06] LocalizeRingForHomalg [BLHIO]. The implementation is abstract in 
the sense that any commutative computable ring R supported by a computer algebra 
system to which the homalg project [hpalO] offers an interface can be localized at any 
of its flnitely generated maximal ideals m, thus providing a new ring R^ for the homalg 
project. The package LocalizeRingForHomalg additionally includes an interface to the 
implementation of MORA's algorithm in SINGULAR [GPS09], which can alternatively be 
used to solve (in)homogeneous linear systems over R^, making (together with the matrix 
operations in the proof of Corollary 4.2) the ABELian category R^ — fpmod computable. 

Our aim to reduce computations over R^ to ones over R suggests the above used well- 
adapted data structure for matrices A over R^: Write A as a fraction ^ with numerator 
matrix A over R and a single denominator s G R\m. Lemma 4.3 and Proposition 4.5, which 
provide the key algorithms for this reduction, require at least common denominators for 
each row of the input matrices. But taking common denominators for each row (or column) 
of a matrix is not suited for matrix multiplication. We saw in the proof of Corollary 4.2 that 
this data structure uses no more than computations of common denominators to realize 
the remaining basic operations for matrices over R^. 

The computational aspects of fraction arithmetic are critical for performance issues since 
writing matrices with a common denominator blows up numerator matrices. It became 
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efficient the moment we started using least common multiples (as far as they exist in R) 
for common denominators and for canceling fractions representing ring elements. 

Solving the submodule membership problem for the same submodule and various dif- 
ferent elements occurs very often in homological computations. For Proposition 4.5 this 



b. Thus replacing A by a distinguished basis (over R, if such a basis exists) is a minor op- 
timization. But with A being a distinguished basis one can ffist check if b reduces to zero 
modulo A. If so, then the row vector z, occurring in the proof of Proposition 4.5, can be 
assumed zero. Hence, the row vector t, occurring in the statement of the proposition, has 
1 as denominator. This heuristic succeeds remarkably often and prevents the creation of 
unnecessary denominators, which would propagate through the remaining computations. 



The examples below are computed using several packages from the homalg project 
[hpalO], all written in GAP4 [GAP06]. Here we use Singular [GPS09] as one of the 
most performant GrOBNER basis engines with an existing interface in the homalg project. 

Example 6.1. Let R := Q[a, b, c, d, e] and M the i?-module given by 4 generators satisfying 
the 9 relations given below, i.e., presented by a 9 x 4-matrix A: 

gap> LoadPackage ( " RingsForHomalg " );; 

gap> R := HomalgFi e IdQf Rat i onal s InSingular ( ) * "a,b,c,d,e"; 
<A homalg external ring residing in the CAS Singular > 
gap> A := HomalgMatrix ( " [\ 

> 2*a+c+d+e-2,2*a+c+d+e-2,2*a+c+d+e-2,0,\ 

> 2*c*d+d~2+2*c*e+2*d*e-3*e~2-2*c-d-6*e+l ,\ 

> 2*c*d+d~2+2*c*e+2*d*e-3*e~2-2*c-d-6*e+l ,\ 

> 2*c*d+d~2+2*c*e+2*d*e-3*e~2-2*c-d-6*e+l ,0,\ 

> -4*a+2*b-c -d-e+2 ,-4*a+2*b-c-d-e+2,\ 

> -c+d+e+2,4*a*d-2*b*d+2*d"2+2*d*e,\ 

> c-2-d-l , c-2-d-l , c-2-d-l ,0 ,\ 

> 4*d*e~2-d"2+2*c*e+4*d*e-3*e~2-2*c+3*d-6*e+5 ,\ 

> 4*d*e"2-d~2+2*c*e+4*d*e-3*e~2-2*c+3*d-6*e+5 ,\ 

> 4*d*e~2-d~2+2*c*e+4*d*e-3*e"2-2*c+3*d-6*e+5,0,\ 

> 0,b"2+a+c+d+e,0,b~2*e+a*e+c*e+d*e+e~2,\ 

> , b-2*d+a*d+c*d+d~2+d*e , , , \ 

> , a*b~2+a~2+a*c+a*d+a*e , , , \ 

> 4*b-3*d-4*d~3-12*d-2*e-32*c*e~2+12*e-3+21*d"2\ 

> -42*c*e+40*d*e+27*e-2+8*c-9*d+16*e-17 ,\ 

> -4*a*b*d-4*b*c*d-4*b*d"2-4*d~3-4*b*d*e-12*d~2*e-32*c*e~2+\ 

> 12*e-3 + 21*d"2 -42* c*e+40*d*e + 27*e-2 + 8*c -9*d + 16*e -17 , \ 

> -12*d-2*e -32* c*e "2+12*6-3 + 21* d~2 -42* c*e+44*d*e+\ 

> 27*e-2+8*c-9*d+16*e-17, -4* b "3+4* d ~2+4* e\ 

> ] " , 9 , 4 , R ) ; 

<A homalg external 9 by 4 matrix> 
gap> M := Lef tPr e sent at ion ( A ); 

<A left module presented by 9 relations for 4 generators> 



means 




many computations with different rows 



6. Examples 
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Let Rq := i?m denote the localized ring at the maximal ideal m = (a, b, c, d, e) in R 
corresponding to the "origin" in A^(Q). We now want to compute (following [Bar]) the 
purity filtration (=equi dimensional filtration) of the localized i?m-niodule 

Mo ■.= M^ = R^®R M. 

One possible approach would be to compute the purity filtration for the global i?-module 
M and to localize the resulting filtration afterwards^^. But since the module M is also 
supported^'^ at components not including the origin it is clear that computing the global 
purity filtration will automatically accumulate any structural complexity of M at these 
components as well. The algorithm suggested in [Bar] for computing the purity filtration 
starts by resolving M. And indeed, an early syzygy computation during the resolution of 
M failed to terminate within a reasonable time. 

The computation for the localized module Mq over the local ring i?o did not terminate 
using Mora's algorithm either, where it gets stuck in a basis computation at an early 
stage. 

However, the purity filtration can be computed for the localized module Mq using the 
approach suggested in §4 within seconds: 

gap> LoadPackage ( " LocalizeRingForHomalg " );; 
gap> RO := LocalizeAtZero ( R ); 
<A homalg local ring> 
gap> MO := RO * M; 

<A left module presented by 9 relations for 4 generators> 
gap> ByASmall erPr e s ent at i on ( MO ); 

<A left module presented by 10 relations for 3 generators> 
gap> filtO := PurityFiltration ( MO ); 

<The ascending purity filtration with degrees [ -3 .. ] and graded parts: 
0: <A zero left module > 
-1: <A cyclic reflexively pure codim 1 left module presented by 

1 relation for a cyclic generator > 

-2: <A reflexively pure codim 2 left module presented by 7 relations for 

2 generators> 

-3: <A cyclic reflexively pure codim 3 left module presented by 

3 relations for a cyclic generator > 
of 

<A non-pure codim 1 left module presented by 10 relations for 3 generators>> 

The following command computes an explicit isomorphism between a new module (equip- 
ped with a triangular presentation compatible with the purity filtration) and our original 
module Mq. 

gap> m := I somorphi smDf Fi It r at i on ( filtO ); 
<An isomorphism of left modules > 
gap> FilteredModule := Source ( m ); 

<A left module presented by 7 relations for 4 generators> 
gap> Display ( FilteredModule ); 
_ [1 ,1] ,0, 0, _ [1 ,4] , 



'Justifying this statement is left to the reader. 
'Recall, supp7\f := {p G Spec(i?) | A/p ^ 0}. 
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0, _ [2,2] ,0, 0, 

0, 0, _[3,3],0, 

0, _[4,2] ,_[4,3] ,-b-l/2, 

, , , _ [5 ,4] , 

, , , _ [6 ,4] , 

0, 0, 0, _[7,4] 

/(4*a*b*d~3-2*b~2*d~3+2*b*d~3*e+2*a*d~3-b*d~3+d~3*e-4*a*b*d+2*b~2*d+2*b*d~2-2*\ 
b*d*e-2*a*d+b*d+d"2-d*e-2*b-l) 

Cokernel of the map 

R"(lx7) --> R~(lx4), ( for R := Q [a , b , c , d , e] _< a, b, c, d, e > ) 
currently represented by the above matrix 



Singular suppressed the relatively big entries of the triangular presentation matrix. We 
use the following command to make them visible. We see that in fact all fractions cancel 
except for one entry which retains a denominator: 



gap> EntriesQf HomalgMatr ix ( MatrixQf Relations ( FilteredModule ) );; 

gap> ListToListList( last, 7, 4 ); 

[ [ (b-2+a+c+d+e)/l , 0/1, 0/1, -1/1 ], 

[ 0/1, (2*a-b+d+e)/l , 0/1, 0/1 ], 

[ 0/1, 0/1, (b-3-d-2-e)/l , 0/1 ], 

[ 0/1, (2*b~2*d*e+b*d*e+2*c*d*e+d"2*e+d*e~2)/l , 

(2*a*b*e+b~2*e+2*b*c*e+2*b*d*e+2*d"2*e+2*b*e~2+a*e+c*e+d*e+3*e~2)/l, 
( -1/2) / (2*a*d-3-b*d~3 + d-3*e -2*a*d+b*d+d-2-d*e -1) ] , 

[ 0/1 , 0/1 , 0/1 , d/1 ] , 

[ 0/1 , 0/1 , 0/1 , a/1 ] , 

[ 0/1, 0/1, 0/1, (b-4-b-3*e-b*e+e-2)/l ] ] 



Remark 6.2. One would wonder why local computations (as in Example 6.1 above and 
Example 6.3 below) using the approach suggested in §4 could be faster than the global ones, 
although the local syzygies computation in Lemma 4.3 is nothing but a global syzygies 
computation and the effective solution of the local submodule membership problem in 
Proposition 4.5 is again nothing but the effective solution of an adapted global one. This 
has to do with the fact that the local ring R^^ has many more units than R, leading to 
the structural simplifications mentioned in Remark 4.4. Furthermore, DecideZerOij^(B, A) 
often equals zero even if DecideZero/j(B, A) (for the corresponding numerator matrices 
B, A) does not (cf. §3.1.2). For the geometric interpretation of the last statement let (C) 
denote the submodule of the free i?-module generated by the rows of the matrix C. Then 
the maximal ideal m often enough does not lie in the support of the nontrivial i?-subfactor 
module (B)/(A), i.e., the i?m-subfactor (B)/(A) is trivial although its global counterpart 
(B)/(A) is not. So one can roughly say that zero modules occur more frequently in local 
homological computations than in global ones, making the former faster, in general. 
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Example 6.3. Serre's intersection multiplicity formula of two ideals J, J<i? at a prime 
ideal p<i? (cf. [Har77, Thm. A.1.1]) 

J; p) = 5^(-l)nength (Xorf^ {R,/I,, RJJ,)) 

i 

offers a nice demonstration of Remark 4.8. 

Let R := ¥^[x, y, z, v, w] with maximal ideal p = xn = {x,y, z, v, w). We use the package 
LocalizeRingForHomalg to define the two localized rings Rq = Sq := R^- The ring 
Rq utilizes Lemma 4.3 and Proposition 4.5, whereas 5*0 uses MORA's algorithm to solve 
(in) homogeneous linear systems. 

gap> LoadPackage ( "Sheaves" ); 

gap> R := HomalgRingOf Integer s InSingular ( 5 ) * " x , y , z , v , w " ; ; 

gap> LoadPackage ( "LocalizeRingForHomalg" ); 

gap> RO := LocalizeAtZero ( R );; 

gap> SO := Local izePolynomialRingAt Zer oWi thMor a ( R );; 



The ideals / and J are the intersection of ideals similar to those in [Har77, Exam- 
ple. A.LLl] with ideals not supported at zero: 



gap> il := HomalgMatrix ( "[ \ 










> x-z , \ 










> y-w \ 










> ] " , 2, 1 , R ) ; ; 










gap> 12 := HomalgMatrix ( "[ \ 










> y~6* v~2*w-y ~3* v*w~20 + l , \ 










> x*y"4*z"4*w-z~5*w~5+x~3*y*z~2-l \ 










> ] " , 2, 1 , R ) ; ; 










gap> I := Intersect ( Lef tSubmodule ( 


il 


) , Lef tSubmodule ( 


12 ) ) 




gap> 10 := RO * I; 










<A torsion-free left ideal given by 


4 


generators > 






gap> 010 := FactorOb j ect ( 10 ); 










<A cyclic left module presented by 


yet 


unknown relations 


for a 


cyclic generator > 


gap> jl := HomalgMatrix ( "[ \ 










> x*z , \ 










> x*w , \ 










> y*z, \ 










> y *w , \ 










> v~2 \ 










> ] " , 5 , 1 , R ) ; ; 










gap> j2 := HomalgMatrix ( "[ \ 










> y~6* v~2*w-y ~3* v*w~2+l , 


\ 








> x*y~4*z"4*w-z~5*w~5+x~3*y*z~2-l , 


\ 








> x-7 


\ 








> ] " , 3 , 1 , R ) ; ; 










gap> J := Intersect ( Lef tSubmodule ( 


jl 


) , Lef tSubmodule ( 


j2 ) ) 




gap> JO := RO * J ; ; 










gap> DJO := FactorOb j ect ( JO ); 










<A cyclic left module presented by 


yet 


unknown relations 


for a 


cyclic generator > 
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Computing the Tor-modules over Sq (using MORA's algorithm) or globally^^ over R both 
did not terminate within a week. But over Rq (using the approach suggested in §4) the 
computation terminates in few seconds: 



gap> TO := 


Tor( 010 


, OJO ); 




<A graded 


homology 


object consisting of 3 left modules at degrees [ . 


. 2 ]> 



As mentioned in Remark §4.8, our approach cannot produce a "distinguished basis" 
for the presentation matrices of the resulting Tor-modules, from which their HiLBERT 
series can be read off. The sum of the coefficients of the HiLBERT series of the module 
Torf'"(i?tn//ni, i?^l/<^a) IS nothing but its dimension as a R/m = F5 vector space, which 
in this case coincides with the length. But now we can apply MORA's algorithm to the 
already computed presentation matrices of the modules Toif"' {Rm/ Im, Rm/ Jm)'- 



gap> TOMora 


:= SO * TO; 








<A sequence 


containing 2 morphisms of 


left modules at 


degrees [ . 


. 2 ]> 


gap> List ( 


Ob j ectsOf Complex ( TOMora 


), AffineDegree 


); 




[ 6, 2, ] 











Thus the intersection multiplicity at m is 6 — 2 + = 4. 

Of course, getting rid of the irrelevant^^ primary components of / and/or J would simplify 
the computations. But we were not able to compute a primary decomposition for / with 
the computer algebra system SINGULAR, and computing one for J took more than seven 
minutes, which is still longer than the few seconds needed by our approach. 



7. Conclusion 

The above axiomatic setup for algorithmic homological algebra in the categories of 
finitely presented modules, merely viewed as ABELian categories, only requires solving 
(in) homogeneous linear systems and does not enforce the introduction of any notion of dis- 
tinguished basis. This abstraction motivated a constructive approach to the homological 
algebra of such module categories over commutative rings localized at maximal ideals, an 
approach in which global computations replace local ones. 

For local polynomial rings this shows that BuCHBERGER's algorithm provides an alter- 
native to Mora's algorithm, an alternative which for some examples may even lead to 
a remarkable gain in computational efficiency. Mora's algorithm remains indispensable 
when it comes to computing HiLBERT series. It is hence a mixture of both algorithms that 
proves more useful in practice. 

Appendix A. Existential Quantifiers in Abelian Categories 

The aim of this appendix is to recall the axioms of an ABELian category. We will view 
them as basic constructions which suffice to build all the remaining ones. In case these 
few basic constructions are computable, it follows that all further constructions become 

^'^Toif'^{S-H'I,S-'^N) = 5-iTorf(M,iV) for multiplicatively closed subsets S C R, cf. [Rot09, 
Prop. 7.17 or Cor. 10.72]. 

"'^^I.e., those primary ideals q with q ^ m. 
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computable as well. Except for the definition of ABELian categories the text closely follows 
[HS97, §11.9]. We will be emphasizing all relevant existential quantifiers which for the sake 
of computability need be to turned into constructive ones. 

We begin by recalling the most basic definitions. So let A he a. category. If not stated 
otherwise we use the following convention^^ for composing morphisms: 

Hom^(M,L) X Hom^(L,A^) ^ Hom^(M,A^), 
{(f) , tfj) ^ 

We call it the left convention, where the right convention would be to write instead 
of (ptp. In any case we call the pre-morphism and the post-morphism. 

A.l. Monies, epics, sums and products. 

Definition A.l ([HS97, §11.3, p. 48]). 

(e) A morphism n : M ^ N is epimorphism (or an epic) if it is pre-cancelable, 
i.e., 

71(f) = Tllf) =^ (j) = l/j 

for all objects L and all morphisms 0, € }lomA{N, L). We write n : M ^ N. 
(m) A morphism l : M ^ N is monomorphism (or a monic) if it is post-cancelable, 
i.e., 

(j)L = IpL =^ (f) = Ip 

for all objects L and all morphisms 0, ^/^ G Homyi(L, M). We write l : M ^ N . 

Definition A.2 ([HS97, §11.5, p. 58, p. 54]). Let / be an index set. 

(p) An object HiG/^* together with a family of epimorphisms {vTj : ~^ is 

called "the" product of if the following universal property is satisfied: 

For any object M and any family {0j : M — ?■ Aj}jg/ of morphisms there exists a 
unique morphism 

called the product morphism, such that ^VTj = 0j for all i E I. 
(s) An object Ai together with a family of monomorphisms {i, : Ai )■ @ Ai}i^i 

is called "the" coproduct of {Ajjjg/ if the following universal property is satisfied: 
For any object M and any family {(j)i : Ai — M}jg/ of morphisms there exists a 
unique morphism 

0=(0,):0A,^M, 
called the coproduct morphism, such that ij0 = 0j for all i E I . 

The universal properties imply that the product and the coproduct, in case they exist, 
are unique up to isomorphism in A. 



16i 



This differs from the convention foUowed in [HS97, §11.1, p. 41]. 
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A. 2. Categories with zeros, kernels, and cokernels. 

Definition A. 3 (Category with zero and zero morphisms [HS97, §11.1, p. 43]). 



consists 



of 



• A zero object is an object in A such that ^1(0, M) and A(M^ 0) 
a single element for all objects M & A (i.e., an object which is both initial and 
terminal). Trivially, such an object is unique up to isomorphism in A. 

• Define for each pair of objects M, N the unique zero morphism 

Omn := M ^ N. 

We write instead of Omn when no confusion is possible. 

Definition A.4 ([HS97, §11.6, p. 61]). Let ^ be a category with and : M A^ a 
morphism. 

(k) A morphism k : K ^ M is called "the" kernel of : M — )■ A^ if 

(i) K,(f) = 0, and 

(ii) for all objects L and all morphisms t : L ^ M with T(f) = there exists a 
unique morphism Tq : L K, such that r = TqK. Tq is called the lift of r 
along K,. 

It follows from the uniqueness of the lift Tq that k is a monomorphism. 




K is called "the" kernel object of Depending on the context kerc/) sometimes 
stands for the morphism k and sometimes for the object K. 
(c) A morphism e : M — )■ C is called "the" cokernel of : M — )■ A^ if 

(i) (pe = 0, and 

(ii) for all objects L and all morphisms t] : N L with (pi] = there exists a 
unique morphism rjo : C ^ L, such that rj = erjo. r/o is called the colift of rj 
along e. 

It follows from the uniqueness of the colift rjo that e is an epimorphism. 




C is called "the" cokernel object of 0. Depending on the context coker some- 
times stands for the morphism e and sometimes for the object C. 



Kernels are cokernels, in case they exit, are unique up to isomorphism. 
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A. 3. Additive and AsELian categories. 

Definition A. 5 ([HS97, §11.9, p. 75]). An additive category ^ is a category satisfying: 



exists 



(Addl) A has a zero object 0. 

(Add2) The producf^ {M®N; tim.t^n) of two objects M and 
(Adds) For all objects M and N in A the set Homyt(M, A^) is an (additively written) 
ABELian group. 

(Add4) For all objects M,N,L in A the composition Homyi(M, A^) x Homyi(A^, L) 
Homyi(M, L) is bilinear. 

Remark A. 6. The following is true for additive categories (cf. [HS97, §11.9, p. 75ff]): 

• The unique zero morphism Omn defined above is exactly the zero element of the 



ABELian group A{M,N). (AddS) requires the existence of an addition and a 
subtraction operation in Homyi(M, A^). 

• For objects M and A^ define the morphisms 

iM ■■= {Im,Omn} ■■ M ^ M ® N and ln := {OnmAn} ■ N M ® N. 

Then ttm'-m + ^'n'^n = lM®Af- 

• It follows that finite coproducts also exist: (M © A^; l^j, t^v) with l^j, Ln as above 
and the coproduct morphism defined by 

(0, ip) := 7rM0 + vttvV^ : M ® N L, 

for two morphisms (p : M L and : N L. 

• In particular, finite products and coproducts "coincide". 

• For K L © M A^ we have {a, ip) = a(j) + Pi). 

• In particular, for (p^ip : M ^ N we have (f) + = {1m, ^n}{<Pi V")- This means that 
the additive structure is determined by the category A. 

In additive categories one often uses the terminology sum instead of coproduct. 

Now we can finally state the definition of an ABELian category: 

Definition A. 7 ([Rot09, §5.5, p. 307]). A category A is called ABELian if 
(Abl) A is additive; 



(Ab2) every morphism has a kernel and a cokernel; 



(Ab3) every monomorphism is a kernel; and every epimorphism is a cokernel. 

It follows that then every monomorphism is the kernel of its cokernel and, dually, every 
epimorphism is the cokernel of its kernel. To see the first statement let k : A' — )■ M be a 
monomorphism and e : M — )■ C its cokernel. To prove that k is the kernel of e, we first 
note that = (as e is the cokernel). Now for a morphism r : L — )■ M with re = 
we need to prove the existence of a unique lift Tq : L ^ K with tqk, = t. To see this let 
1] : M ^ N he a morphism of which k, is the kernel (Ab3). Since e is the cokernel of k 



17i 



The use of the coproduct symbol wiU be justified below. 
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there exists a (unique) colift tjq with erjQ = rj. Hence rrj = rerjQ = Orjo = as well and there 
exists a unique lift Tq : L ^ K such that TqK = r and we are done. 



It also follows that in ABELian categories every morphism is expressible as the com- 
position of an epimorphism (namely vr := coker(ker0)) and a monomorphism^^ (namely 
L := ker(coker 0)). More precisely, the homomorphism theorem in ABELian categories 
takes the following form: 



Proposition A.8 ([HS97, Prop. II.9.6],[ML95, Thm. IX.2.1]). For (p : M ^ N in an 

ABELian category there is a sequence 



where ttl = (j), k is the kernel of (f) and vr, e the cokernel of (p and l, vr is the cokernel of k, 
and L is the kernel of e. 

A. 4. Derived functors and projective resolutions. Functors and their natural trans- 
formations play a central role in category theory. Functors between ABELian (or additive) 
categories preserving direct sums of objects are called additive. A powerful tool in homo- 
logical algebra is the construction of the higher derived functors of a non-exact additive 
functor F-.A-^'B [HS97, Chap. IV]. 

To define the higher left derived functors of a (covariant) functor F one usually requires 
the category to have enough projectives (see Def. A. 10). This is needed to construct projec- 
tive resolutions (see below). For right derived functors replace "projective" by "injective". 
The derivation of a contravariant additive functor G : — )■ 3 is defined by viewing G 
as a covariant functor G : — ?■ 23. A projective resolution in 71°^ becomes an injective 
resolution in A and vice versa (for details cf. [HS97, §11.4, §11.5]). 

In fact there is a definition of derived functors which works even if the source ABELian 
category does not have enough projectives: Let F : A ^ Ab be a right exact (and hence 
additive) functor from the ABELian category A into ABELian groups. Define {LqF)A : = 
[Ext'^(A, — ), F], the ABELian group of natural transformations from the functor Ext''(A, — ) 
to the functor F (cf. [HS97, Cor. IV.10.2]). The point is, that now Ext''(A,E), the value 
of the derived bifunctor Exf, can be described as the ABELian group of YONEDA's q- 
extensions of two objects A,B ^ A without referring to projective or injective resolutions 
(cf. [HS97, Ex. 111.2.5,2.7, §IV.9]). 

There are two drawbacks of this definition. First, the functor F must have values in the 
category Ab of ABELian groups. Second, it is a priori not clear how to compute Ext'^(y4, B) 




L 



VO 



'This is stated as the last defining axiom of an ABELian category in [HS97, §n.9, p. 78]. 
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as the AsELian group of Yoneda g-extensions, let alone the "set" [Ext''(y4, — ), F]. Since in 
this paper we will only be using module categories as an example to illustrate our axiomatic 
approach to computability, we will define derived functors via projective resolutions. More 
precisely, starting from the next section we will focus on the computability in the category 
A := R — fpmod of finitely presented modules over some ring R. It turns out that 
for many rings of interest all injective objects in i? — fpmod are zero objects. In contrast, 
R — fpmod always has enough projectives. 

For this reason we will restrict our discussion to AsELian categories having enough 
projectives. This restriction only allows us to left derive covariant additive functors and 
right derive contravariant ones. 

In rest of this subsection we will recall especially those notions needed to emphasize two 
further existential quantifiers: 



exists 



Definition A. 9. An object P in a category A is called projective, if for each epimorphism 
e : M ^ N and each morphism (f) : P ^ N there 
0ie = 0. 

P 



a morphism 0i : P — )■ M with 



y 

M 




N 



We call 01 a projective lift of cj) along e. 



A supposedly more general form of the projective lift is often used in ABELian categories. 
The assumption of e being epic can be relaxed in the following way: Let A be an ABELian 

category and e : M — )■ a morphism with e : M ^ I N. According to Prop. A. 8 
there exists an essentially unique decomposition of e into an epic vr and a monic l. Further 

let /3 : — 7- L be a morphism with kernel l (in other words, M A^ A L is an exact 
sequence), P & A projective object, and : P — )■ A^ a morphism with 0/3 = 0. This last 
condition expresses that the image subobject^^ of is "contained" in the image subobject 
of e. It easily follows that there exists a projective lift 0i along e making the following 
diagram commutative: 

P ^ 





The projective lift 




is constructed in two steps: 
P 




19 



For details see [Rot09, Def. of suboject, p. 306]. 
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First construct 0o as the lift of • 
of along the epimorphism vr. 



along the monomorphism l. 0i is then the projective lift 



Definition A. 10. A category A is said to have enough projectives, if for each object 
M there exists an epimorphism u : P ^ M with P projective. The epimorphism u is 
called a projective hull of M. The kernel object Ki := kei u of its kernel fii : Ki ^ Pi 
is called the 1-st syzygy object of M. 

Having enough projectives we can again find a projective hull Ui : Pi ^ Ki of the 
first syzygy object Ki. Iterating this process yields the higher syzygy objects Kn, 
embeddings yU„ : t- Pn-i? projective hulls '■ Pn ^ Kn, and a projective resolution 
P, of M with Pq := P and (9„ := Unf^n'- 




M Ki K2 

The morphism 5i : Pi — )■ Pq is called a projective presentation of M = coker(9i. 

Appendix B. Further constructions 

B.l. Further constructions in ABELian categories. There are several others categor- 
ical constructions which automatically exist in AsELian categories. Hence, for establishing 
their existence we only need the basic ones listed in §2. 

We will demonstrate this for the pull-back and (its dual) the push-out, being the two 
most prominent ones. 

Definition B.l. Let yi be a category. A pull-back of two morphisms </> : L — )■ M and 

■j/) : J — )■ M is an object X in A and a pair of morphisms A : X — )■ L and j : X ^ J such 
that X(f) = jip satisfying the following universal property: 

For each object Z and each pair of morphisms p : Z ^ L and f] : Z ^ J with p0 = rjijj, 
there exists a unique ( : Z ^ X with p = (X and 77 = (j. 

(2) 




Reversing all arrows in the above definition we obtain the dual definition of the push-out 
of two morphisms : M — )■ L and ip : M J. 

Next we recall why finite pull-backs and push-outs exist in ABELian categories and how 
they can be constructed using kernels and cokernels, respectively. 
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Remark B.2. In an ABELian category A the following two dual statements hold. 

(1) The pull-back oiL ^ M ^ J oi two morphisms can be computed as a kernel. More 

(<t>,'<P) 

precisely X ^ © J is the kernel of the coproduct morphism L © J . 

The unique morphism C, in diagram (2) is nothing but the lift of {p, — r/} along the 
monomorphism {A, — j}. 

(2) The push-out oi L M ^ J oi two morphisms can be computed as a co- 
kernel. More precisely L ® J — is the cokernel of the product morphism 

m—-Il®J. 

In fact, the existence of finite pull-backs and push-outs is a consequence of the existence 
of finite limits and colimits in additive categories with kernels and cokernels (cf. [Rot09, 
Ex. 5.60, p. 321]). 

B.2. Further constructions in module categories. Categories of modules and cat- 
egories of sheaves of ABELian groups are the most prominent ABELian categories. For 
the connection between ABELian categories and module categories see "the full embedding 
theorem" in [Fre64, Chap. 7] and the discussion in [Har77, §111.1, p. 203]. 

As objects in module categories are sets it is natural to ask whether one can describe 
certain set theoretic operations for modules in a purely categorical way. We will give two 
examples: 

The set-theoretic intersection of two submodules L of an i?-module M can be de- 
scribed as the pull-back^° oi K ^ M L, where lk '■ K ^ M and ll '■ L ^ M are the 
corresponding embeddings. 

The submodule quotient 

K : J := {r e R\rJ C K} 

of two submodules K, J of a left module M is a yet more specific construction in module 
categories. The computation of the left ideal K : J can be expressed in terms of an 
intersection of kernels of certain morphisms 

K -.J = p|ker0^, 

where 0m '■ R ^1 K, r rm for an m G M. Of course, the intersection might be taken 
over a generating set of J. 
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